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ABSTRACT 

Using the first 25% of DEEP2 Redshift Survey data, we probe the line-of-sight velocity dispersion 
profile for isolated galaxies with absolute B-band magnitude —22 < Mb ~ 51og(/i) < —21 at z=Q.7- 
1.0, using satellite galaxies as luminous tracers of the underlying velocity distribution. Measuring the 
velocity dispersion beyond a galactocentric radius of ^ 200ft-~"'^kpc (physical) permits us to determine 
the total mass, including dark matter, around these bright galaxies. Tests with mock catalogs based 
on N-body simulations indicate that this mass measurement method is robust to selection effects. 



We find a line-of-sight velocity dispersion {(Tios) of 1(!>2'^^ km s ^ at ~ 110ft. ^ kpc, I36I20 ^ 



5+26 



at 



230/1^1 kpc, and ISOtgg km s^^ at - 320/1"^ kpc. Assuming an NEW model for the dark 



matter density profile, this corresponds to a mass within r2oo of A/200 ~ 5.5 



+2.5 



IQi^ft-iMm for 



our sample of satellite hosts with mean luminosity ^ 2.5L*. Roughly ^ 60% of these host galaxies 
have early-type spectra and are red in restframe [U — B) color, consistent with the overall DEEP2 
sample in the same luminosity and redshift range. The halo mass determined for DEEP2 host galaxies 
is consistent with that measured in the Sloan Digital Sky Survey for host galaxies within a similar 
luminosity range relative to M^. This comparison is insensitive to the assumed halo mass profile, and 
implies an increase in the dynamical mass-to-light ratio {M2qo/Lb) of isolated galaxies which host 
satellites by a factor of ^ 2.5 from z '--^ 1 to z ~ 0. Our results can be used to constrain the halo 
occupation distribution and the conditional luminosity function used to populate dark matter halos 
with galaxies. In particular, our results are consistent with scenarios in which galaxies populate dark 
matter halos similarly from z ~ to z ~ 1, except for ~ 1 magnitude of evolution in the luminosity 
of all galaxies. With the full DEEP2 sample, it will be possible to extend this analysis to multiple 
luminosity or color bins. 

Subject headings: galaxies: evolution — galaxies: kinematics and dynamics — galaxies: halos — dark 
matter 



1. INTRODUCTION 

It has been firmly established that galaxies and clusters 
form within halos whose mass is dominated by unseen 
dark matter. Yet until very recently, the outer regions of 
halos have been very poorly understood due to a lack of 
visible tracers of the mass distribution. Galaxy-galaxy 
lensing is able to probe the halo masses of local galaxies, 
though with some degree of uncertainty, as this method 
actually probes all of the mass along the line-of-sight 
(jGuzik fc SeliakI 12002(1 . and ha s only recently been ap- 
phed t o isolated gala xies (Hocks tra et alJl2005l). Recent 
work ( Wilson et al. 2001: Hockst ra. Yee. fc Gladder^ 
1200 4: Kle inhcinrich et al . 2005) suggests that the virial 
mass of ~ L* galaxies has remained constant from z ~ 
0.8 to z ~ 0.15. Beyond z ~ 0.5 the lensing probabil- 
ity rapidly diminishes, making it very difficult to derive 
masses of isolated galaxies (with masses ^ lO^'^ Mfr,) with 
this method at higher redshift (see e.g.. iPeacock 1999). 

The dynamics of satellite galaxies orbit- 
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ing larger "host" galaxies provide another 
way to probe the ma ss distri bution at large 
radii. Early work bv ILittle fc Trcmainc ( 19872; 
'Erickson. Gott esman. fc Hunted l|1987i|l : iZaritskv et al] 
(1989) utihzed samples of tens of satellites as early 
confirmation that galaxies are embedded in large dark 
matt er halos. M ore conclusive evidence was compiled 
by Z aritskv et al.l (^993, 199^ who used the kinematics 
of a sample of 115 satellite galaxies to probe the outer 
regions of 69 isolated galaxies. By employing satellites as 
test particles, they built up a velocity profile for a single 
representative isolated galaxy by stacking measurements 
of satellites from many different hos t galaxies. 

More recently, 'Prad a et all l|2003[ hereafter P03) use 
^250, 000 SDSS redshifts to probe the halo masses of iso- 
lated galaxies; they detect >1000 satellites around '~700 
host galaxies. With this large data set they are able to 
discriminate between various halo mass distributions and 
find evidence for an NFW-like falloff {p oc r~^) at large 
radii. From these accurate line-of-sight velocity disper- 
sion profile measurements, P03 infer the masses enclosed 
within 1.5 Hviriai for two sets of isolated galaxies. Host 
galaxies with —20.5 < Mb < —19.5 are found to have 
an average halo mass of Mmriai ~ 1-5 x 1O^^M0 while 
hosts with —21.5 < Mb < —20.5 have Myiriai ~ 6 x 
lO^^M© (for h = 0.7). Other recent work u tilizing satel- 
lite dynamics includes iMcKav et al.l l)2002ii) . who check 



2 



SDSSweak-lensing scaling laws; Ivan den Bosch et al.1 
l)20n4alhD . who use mock galaxy catalogs and the 2dF 
survey to constrain the conditional luminosity function 
and investigate the levels and effect s of contami nation 
in dynamical satellite studies; and Braincrd (2005), who 
measure velocity dispersion profiles for subsamples of iso- 
lated 2dF galaxies. 

By extending this type of measurement to high red- 
shift, we can study the evolution of the relationship 
between galaxies and dark matter halos. There have 
been few ways to do this until the recent advent of 
large, high-redshift survey s . The best example to date is 
lYan. Mad gwick. fc Whit3 l|2003j) . who use the 2dF and 
DEEP2 t wo-point correlation functions l|Madgwick et al.| 
l2003at iCoil et all 12004 respectively) to constrain the 
evolution of the halo occupation distribution, a key in- 
gredient of the halo model. Yan et al. produced a set of 
z ~ 1 mock catalogs using N-body simulations and a halo 
model whose parameters are set by requiring a fit to ^(r) 
from 2dF at z ~ 0. They find good agreement between 
the correlation statistics at z ~ 0.8 from DEEP2 and 
the prediction from these mock catalogs, which populate 
galaxies in dark matter halos in the same way (as a func- 
tion of L/L* and halo mass) at z 1 and z ~ 0. Hence 
their results are consistent with a minimal-evolution hy- 
pothesis, in which galaxies with a given luminosity com- 
pared to L* at z ^ 1 live in the same sorts of halos as 
similar galaxies at z ^ 0, though the mass function of 
dark matter halos and L* evolve. Here we address this 
hypothesis with an independent method. 

In this paper we constrain the velocity dispersion pro- 
file for a typical isolated DEEP2 galaxy at z ~ 0.8 us- 
ing methods similar to those of P03. We then deduce a 
representative halo mass for these galaxies and compare 
our results to recent local measurements from SDSS to 
test for evolution in a self-consistent way. We use mock 
catalogs to test the significance and robustness of these 
results. The paper proceeds as follows. In § |21 we de- 
scribe the DEEP2 data set and the properties of satellite 
galaxies and their hosts. § |31 outlines our method for 
reconstructing the mass of isolated galaxies using satel- 
lites and in §^we present our results and compare with 
recent local measurements. We test our methods using 
mock catalogs in § [3 and discuss some implications of 
our results in § El Throughout the paper we assume a 
standard ACDM cosmology with = 0.3, JIa = 0.7 
and Ho = 100/1"^ km s^^ Mpc'^ Absolute magni- 
tudes quoted have been K-corrected and corrected for 
reddening by galacti c dust, and are in the AB system 
^Willmer et ablboO,^. 

2. SATELLITE-HOST SYSTEMS AT Z ~ 1 

In this Section we introduce the data used at z ^ 1, 
describe the algorithm used to identify bright isolated 
galaxies and their satellites, and highlight several prop- 
erties of these host-satellite systems. 

We use data from the first ~ 25% of the DEEP2 
Galaxy Redshift Survey, a three-year project using the 
DEIMOS spectrograph at the 10-m Keck II telescope to 
survey galaxies at z ~ 1. DEEP2 will collect spectra of 
^50,000 galaxies from 0.7 < z < 1.4 to a limiting mag- 
nitude of Rab = 24.1 with redshift errors of ^ ^20 km 
. For survey details, see iDavis et al.l l)2004j) . Pho- 
tometric data were taken in the R and / bands with 



the CFH12k camera on the 3.6-m Canada-France-Hawaii 
telescope. We use data taken during the first two sea- 
sons of DEEP2 observations, which have yielded ~12,000 
secure redshifts over ~0.9 sq. degrees. Our observed R- 
band limiting magnitude corresponds to a different rest- 
frame wavelength with redshift, from 4000 A at z = 0.7 
to 2800A at z = 1.4. This results in a different selection 
function for red and blue galaxies with redshift, such that 
as we move to fainter magnitudes, red galaxies become 
undetectable before blue galax ies; this effect in creases 
with increasing redshift (see Willmer et alJl2005l for de- 
tails). To minimize this effect we consider only bright 
hosts with z < 1. 

We define an isolated galaxy as having no bright neigh- 
bors within a given search cylinder. Once isolated galax- 
ies have been identified, we use another search cylinder to 
identify faint satellite companions. Specifically, a galaxy 
is isolated if it has no neighbors in the DEEP2 spec- 
troscopic sample within a physical distance projected on 
the sky Vp < 350/i~^kpc, linc-of-sight velocity difference 
|Aw| < 1000 km s~^ and absolute magnitude difference 
AMb < 1.5. An isolated galaxy furthermore cannot 
have any neighbors within 350ft.~^kpc< rp < 700/i~^kpc 
and \Av\ < 1000 km s"^ with AMb < 0.75; we relax 
our magnitude cut at large rp because galaxies this far 
apart will be less dynamically associated. Satellites are 
similarly defined to be galaxies with rp < 350/i^^ kpc, 
\Sv\ < 500 km s~^ and 6Mb > 1-5; i.e. satellites must be 
more than 1.5 magnitudes fainter than the host galaxy 
they belong to. These parameters define sample 1 in 
Table 1 which lists the search parameters used in this 
analysis along with the number of found satellites and 
several derived host halo parameters. We consider 7 dif- 
ferent search criteria and find that the results presented 
are quite insensitive to variations in these criteria; for 
convenience we quote results from sample 1 unless oth- 
erwise noted. Our choices of parameters ensure that a 
satellite can be associated with one and only one host 
galaxy. We put no restriction on morphological or spec- 
tral type, but require the isolated host galaxy to have 
-22 < Mb - 51og(/i) < -21 and 0.7 < z < 1.0. 

For our chosen set of search parameters (Sample 1), 
we have identified 75 satellites around a total of 61 host 
galaxies at 0.7 < z < 1.0 in the DEEP2 data set. Fig- 
ures n and [3 show relevant characteristics of the satel- 
lite and host galaxies including distributions in redshift, 
spectral-type, absolute magnitude, satellite number per 
host, satellite distance from host, and AMb between the 
host and satellite. We determine spectral types using the 
principal component analysis of Madgwic k et alJ l)2003bf) 
and use their definition ofrj = —13 to separate early and 
late-type galaxies. Galaxy morphology and {U ~ B)o 
color correlate well with this spectral classification of 
early and late types l^adgwick 2003)). Satellites are 
found to have ^ 90% late-type spectra, but due to the 
DEEP2 survey selection effects mentioned above, it is 
difficult to determine if this is an intrinsic property of 
satellites around bright isolated galaxies, or due to the 
i?-band selection of the survey. As we observe fainter 
galaxies (e.g. satellites), early- type galaxies become un- 
detectable before late-type galaxies. 

We find that approximately 60% of host galaxies have 
early-type spectra. When we select a subsample of the 
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Fig. 1. — Characteristics of isolated host galaxies. Upper left: 
absolute B-band magnitude of hosts. Upper right: redshift distrib- 
tion. Lower left: first PCA eigenvalue for host galaxy spectra. The 
dashed line indicates the division between late types {r) > —13) and 
early types (77 < —13). Lower right: number of satellite galaxies 
found per host. 
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Fig. 2. — Characteristics of satellite galaxy properties. Upper 
left: satellite galaxy absolute B-band magnitude. Upper right: 
AMs between host and satellite. Lower left: satellite galaxy spec- 
tra type. The dashed line indicates the division between late and 
early type spectra. Lower right: projected separation (physical 
units) on the sky between satellite and host. Slit collisions reduce 
the number of satellites found with Tp < 50/i~^ kpc. 



entire available DEEP2 data set such that it has the 
same redshift and absolute magnitude distributions as 
the host galaxies, we find that both sets of objects 
have consistent fractions of early- type galaxies: 58% 
for isolated hosts and 63% for the subsample. This 
result is somewhat surprising; one might have naively 
suspected that the majority of early-type galaxies with 
—22 < Mb — h\og{h) < —21 would reside in dense envi- 
ronments and hence would not be identified as isolated 
using our search criteria. To test the robustness of this 
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Fig. 3. — Color-Magnitude diagram for all isolated host galaxies 
from sample 1 (stars) compared to a subsample of the first two 
seasons of DEEP2 data (points). The subsample was constructed 
to have the same redshift distribution as the host galaxies, thus 
reducing any potential redshift dependent selection effects. Dashed 
lines indicate the magnitude range we have used in computing halo 
masses. We classify galaxies with (U — B)o >1 "red" and (U — 
B)o <1 "blue". 



result, we also use (U — B)q color to test for any dif- 
ferences between isolated galaxies and our reconstructed 
subsample. Again we find that ~ 55% of isolated host 
galaxies are red {{U — B)o >1), while ~ 50% of the sub- 
sample is red (see Fig. There thus seems to be a 
significant population of bright, red, early- type isolated 
galaxies with satellites at z 1. Our host galaxies have 
approximately the same relative number of early and late 
spectral types as in P03's sample; this will allow for ro- 
bust comparisons between halo masses using SDSS data 
at z - and DEEP2 data at z 1. 

3. HALO MASS ESTIMATION 

In this section we describe our method for obtain- 
ing dark matter halo masses for isolated galaxies uti- 
lizing satellite galaxy kinematics. Our approach is 
similar to previous work (see e.g. IPradaet all 12001 
iBrainerd & Soecian 200^ . except for our maximum- 
likelihood approach to deriving velocity dispersions 
which is more robust than previous methods when ap- 
plied to small numbers of host-satellite pairs. 

Schematically, we derive masses in the following way. 
First we obtain a sample of isolated galaxies with as- 
sociated satellites. We then measure the line-of-sight 
velocity dispersion in several radial bins for a "typical" 
isolated bright galaxy, using satellites as luminous trac- 
ers of the velocity field. In our sample, each isolated 
galaxy has at most three or four satellites, but we can 
measure dispersions by stacking the satellite-host pairs 
and thereby treating all satellites as belonging to a single 
typical host galaxy. In implementing this method we are 
assuming that similarly bright isolated galaxies reside in 
similar halos. We create an homogeneous host galaxy 
sample by searching for satellites around galaxies with 
-22 < Mb - 51og(ft,) < -21 and 0.7 < z < 1.0. As 
more data becomes available, it will be possible to limit 
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these criteria even further, yielding results for multiple 
luminosity, redshift, color, and spectral type bins. 

In order to determine line-of-sight velocity dispersions, 
we build distributions of the projected velocity difference 
between satellite and host (6v) in bins of projected ra- 
dius from the host galaxy and fit for the dispersion in the 
distribution. To convert the resulting velocity dispersion 
measurements into a mass, we make several assumptions, 
including the shape of the underlying dark matter poten- 
tial. We now describe this procedure in detail. 

3.1. Velocity Dispersion Measurement 

The difference between the host and satellite galaxy 
line-of-sight velocity, Sv, has a distribution that is well 
fit by a Gaussian of zero mean. In the absence of inter- 
lopers (see below), the satellite velocity dispersion, aios 
can simply be obtained from the width of a Gaussian 
fit to the velocity distribution. In order to be able to 
detect variation in the velocity profile with radius, we 
bin the satellites in projected radius, rp, from the host 
galaxy. We have chosen bins such that the number of 
satellites per bin is roughly constant: 30 < < 180, 
180 < Tp < 280, and 280 < rp < 350, in units of 
kpc; these choices assure similar errors from bin to bin. 

An important aspect of this analysis is the careful re- 
jection of "interlopers" ; these are galaxies which meet 
the criteria for satellite identification, but are in fact not 
dynamically associated with the host. Interlopers result 
from peculiar velocities which can, in redshift space, scat- 
ter objects into our search cylinder. Re cent local studies 
(e.g. P03. Ivan den Bosch et alJl20n4ajl have found that 
about 20-30% of putative satellites fall into this category. 

Since interlopers are not physically associated with the 
host galaxy, we account for them by assuming that they 
will have a constant Sv distribution. Thus, we expect 
the observed 6v distribution be a combination of flat 
(interloper) and Gaussian (true satellite) components. 
We therefore fit a Gaussian plus constant distribution 
to the velocity measurements within each rp bin. If we 
ignore clustering effects, which is reasonable since we are 
only probing isolated systems, then one would expect the 
number density of interlopers to simply scale with the 
search volume. Although the interloper fraction should 
be roughly constant in Sv, that constant should be differ- 
ent in bins of different radii. This is another important 
motivation for measuring aios in bins of rp. 

Unlike previous studies which fit Gaussian profiles 
to velocity histograms, we determine the dispersion of 
the velocity distribution using a maximum likelihood 
Gaussian-plus-constant fit to the unbinned 6v data. 
Specifically, our likelihood function: 



(1) 



has two free parameters, a constant component (in satel- 
lites per km/s), a, and the width of the Gaussian, aios- 
The parameter B is chosen such that the integral of the 
probability density function of the relative velocity dis- 
tribution between host and satellite (given a and aios) 
over the allowed velocity range is one, and Svi is the Sv 
for the ith satellite-host pair. We maximize the summed 
logarithm of this likelihood function: 



S{a, aio 



Y,ln{L) 



(2) 



over a dense grid in aios and a. 

In Monte Carlo tests, this algorithm agrees with 
method of fitting Gaussian distributions to velocity his- 
tograms for well-sampled data, but is much more robust 
in the limit of small numbers of satellites. This technique 
also provides an estimate of the error in the velocity dis- 
persion measurement from the width of the likelihood 
peak. 

3.2. Halo Mass Determination 

In order to derive a host halo mass, we fit a theoret- 
ical velocity dispersion profile to the measured velocity 
dispersion points. The theoretical profile is obtained via 
the following procedure. We sta r t by a ssuming an NFW 
l|Navarro. Frenk. fc Whit^llQQfl 11997!) density distribu- 
tion 



p(r) 



Sc 



(3) 



{r/rs){l + r/rsf 

(p oc r~^ for large r) where is the present critical 
density, = r2oo / c, and 



Sr.= 



200 



3 ;n(l + c) -c/(l-f c) 



(4) 



where r2oo , defined as the radius where the mean interior 
density is 200 times the critical density. The concentra- 
tion, c, can be viewed as a free parameter determining 
the shape of the NFW density profile. In general the 
concentration is inversely related to the mass of a dark 
matter halo. 

The Jeans equation is then used to relate the radial 
velocity dispersion, ct^, to the gravitational potential, 



- — {pal) + 213^ = - — 
p ar r ar 

and then we integrate along the line of sight 



f^fosirp) 



where 



^M{rp) = 2 



(1-/3-1) 



rp{r 



rl^paf{r,f3)r 



(5) 



dr (6) 



=dr 



(7) 



is the surface mass density fsee lLokas Sz MamonI ll2001fl 
for details of these calculations). In the above, r is the ra- 
dial distance and Vp IS, 0(S usual, the distance projected 
on the sky. The velocity anisotropy (/3 = 1 — al/a'j_, 
where ar is the radial velocity dispersion and a± is the 
velocity dispersion perpendicular to the line of sight) 
must be assumed in the conversion of the NFW den- 
sity profile to a velocity dispersion profile; we use an 
Osipkov-Merrit anisotropy , Pnu — s^/js"^ + s i ), with 
Sa = 4/3 and s = r/ro nn. Ivan den Bosch et alJ l|2nn4a|) 
and (iMamon &: Lokaal2005D find that the line-of-sight 
velocity dispersion profile depends only weakly on (3, at 
the level of a few percent, and hence we do not explore 
other parameterizations. We further need to assume the 
concentration; we use c = 10, which is consistent with 
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our fit to the mock catalog velocity dispersion profiles 
(see §inil, and is in general in agreement with simulations 
of ^ IO^^Mq halos. In §0]we show that our results con- 
cerning the evolution of the halo mass of isolated galax- 
ies are insensitive to the assumed concentration. With 
these assumptions we are left with only one free param- 
eter, the normalization of the velocity dispersion profile, 
which can be characterized by the circular velocity at 
'''200 J ^^200- We fit for V200 via minimization using the 
observed data points. For the same enclosed region, the 
interior mass (M200) ca n be easily inferred from Vmo for 
a giv en cosmology (see iNavarro. Frenk. fc Whitel ll99fiL 
119971 for details). 

4. RESULTS 

Using the maximum likelihood method outlined in §|3 
we measure a velocity dispersion of 162j]3Q km s^^ for 
satellites with 30 < < 180/i~^kpc (median rp = 
llO/i^^kpc), 136tlo km s'^ for 180 < Vp < 280/i-ikpc 
(median Vp — 230/i^^kpc), and ISOlgg km s^^ for 
280 < Tp < 350/1" ^kpc (median Vp = 320/1" ^kpc) for 
isolated galaxies with —22 < Mb — 51og(/i) < —21 (see 
Fig. Errors on the velocity dispersion are derived 
from the maximum- likelihood fit. These results are ro- 
bust to changes in the search parameters and radial bin- 
ning; for the 7 different search criteria listed in Table 1, 
our line-of-sight velocity dispersion measurements vary 
within IcT of the dispersions quoted above. From Fig. ^ 
it is clear that the derived dispersion profile is consistent 
with nearly all popular ha lo mass density pr ofiles (e.g. 
isothermal, NFW, Moore ijMoore et al.lll998|) ). though 
we use an NFW profile to derive masses. What is impor- 
tant for this analysis is the normalization of the velocity 
dispersion profile, not the slope, as long as the slope at 
V200 is shallow. 

As outlined in §|31 we fit velocity dispersion profiles de- 
rived from an NFW model to the measured velocity dis- 
persion profile for DEEP2 by minimizing ^ . Our results 
imply a total mass, M200, of 5.5j;2;o ^ Vd^'^hr^M^ for a 
typical isolated galaxy with —22 < Mb — 51og(/i) < —21 
and at least one satellite. When we vary the search crite- 
ria the measured mass varies by ±1 X IQ^'^h-^MQ, within 
our la errors (see Table 1). As mentioned above, inter- 
lopers play a key role in this analysis. From our best 
fitting Gaussian plus constant fit to the satellite 5v distri- 
bution, we find that the interloper fraction increases with 
increasing rp, from ^ 7% at Vp ^ llO/i"^ kpc to ~ 30% at 
Tp = 320/1"^ kpc. These numbers a re in good agreement 
with pr evious results from P03 and Ivan den Bosch et al.l 
ll2004aD . 

We can compare this derived mass to recent local mea- 
surements to measure evolution in the halo mass of iso- 
lated galaxies. Unfortunately, we could not implement 
the exact same search criteria as in P03 as that yields 
only 30 satellites in the DEEP2 sample, too few to pro- 
vide robust results. Thus, in order to make a fruitful 
comparison, we have taken the raw A?; and rp measure- 
ments of P03 (F. Prada 2004, private communication) 
and independently determined the underlying halo mass 
using our own host galaxy magnitude range and search 
criteria. We note that our methods accurately recover 
the masses inferred in P03's samples (see §nj when us- 
ing their definitions and absolute magnitude intervals. 
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Fig. 4. — Line-of-sight velocity dispersion profile for a typical 
isolated bright galaxy in the DEEP2 sample. The line with error 
bars is the profile derived from the DEEP2 sample; the dotted 
curve is a best fitting NFW line-of-sight velocity dispersion profile 
(c = 10 assumed). The line with errors is very stable over a wide 
range of search parameters and variation in the radial binning. The 
individual points are the At; and rp of the satellite-host systems. 
The best fit NFW profile (dashed line) corresponds to a mass of 
M200 = 5.5 X lOi2/i-iM0. 
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Fig. 5.— As Fig. ^ but for SDSS satellite galaxies from P03 
(data provided by F. Prada). The best fitting NFW line-of-sight 
velocity dispersion profile implies an average halo mass M200 of 
5.4ti;o X 1O^2/i"1M0. The lack of points past rp = 350/i"^kpc is 
a result of the definition of satellites in the P03 dataset. 



Taking MU z ~ 0) - 51og(/i) -19.45 ± 0.07 
ijNorberg et al.ll2002L this value was converted from the 
6j-band by assuming a median color of B — V = 0.21) 
and M Uz ~ 1) - 5 \og{h) = -20.6 ± 0.1 l|Wi1lmer et alJ 
|2005|) . we use 475 host galaxies in the P03 sample with 
magnitudes —21 < Mb — 51og(/i) < —20 ensuring that 
the samples at 2; ^ 1 and z ~ have similar host galaxy 
magnitude ranges relative to Mg . We find isolated galax- 
ies at z ~ with —21 < Mb — 5 log(/i) < —20 to have an 
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average halo mass of 5.4t|^ Q x lO^^h ^Af© when assum- 
ing an NFW mass density distribution with c = 10 (see 

Fig. 13. 

We can cast these results in terms of the dynamical 
mass to B-band light ratio, M2oo/Lb- Using the mean 
luminosity of host galaxies at z ^ 1 of Mb — 5 log(/i) — 
-21.5 ±0.1 {Lb = 5.3 ±0.5 x lO^^h'^LQ) and at z - of 
MB-51og(/i) = -20.5±0.1 {Lb ^ 2.1±0.1x 10^%"^ Lq), 
we find that the S-band mass-to-light ratio {M2oo/Lb) is 
increasing from M200/LB = 104±43 HMq/Lq^b at z ~ 1 
to M200/LB = 257 ± 54 HMq/Lq^b at z - 0, a factor 
of 2.5. Total (random and systematic) errors in our ab- 
solute magnitude measurements, including uncertainties 
in K corrections apart from the assumed cosmological 
parameters, are esti mated to be below .1 mag at the 
redshifts of interest l|Willmer et al.ll2005D : hence errors 
in luminosities are negligible compared to the statistical 
uncertainties in the average M2oo/Lb- 

5. TESTING THE METHOD 

We use mock galaxy catalogs that have been con- 
structed to match the DEEP2 survey in order to test 
whether we can accurately recover the halo mass of iso- 
lated galaxies using satell ites. A complete description of 
the catalogs are given in lYan et all ^^2004^ ; we give the 
relevant details here. N-body simulations of 512^ dark 
matter particles with a particle mass rupart = 1.0 x 10^" 
h~^M© were run in a ACDM Universe using the TreePM 
code (White 2002) in a periodic, cubical box of side 
length 256 h^^ Mpc. Dark matter halos were identified 
by running a "friends-of-friends" halo finder with a min- 
imum size of 8 dark matter particles. Galaxies down to 
O.IL* are then assigned to individual dark matter parti- 
cles via a halo model approach, in which both the number 
of galaxies within a halo and th eir luminosity function 
depends upon halo mass l|Yang et al. 2003). 

We investigate the effects of slitmask target selection 
(see lDavis et al.l2004 for details) on the derived velocity 
profile, thus testing our ability to recover the true halo 
mass when these observational effects are included. Tar- 
get selection may result in a galaxy being identified as 
isolated which is not truly so, since we only have red- 
shifts for approximately one-half of the galaxies meet- 
ing the survey selection criteria. Eventually, photometric 
redshift estimates of galaxies without spectroscopy will 
allow us to better constrain the number of truly isolated 
galaxies in our sample. Earlier work (P03) has inves- 
tigated interloper contamination using simulations of a 
single IO^^Mq halo; the large cosmological simulations 
we use also test the effects of higher order clustering on 
the interloper problem. 

Fig. El compares a line-of-sight velocity dispersion pro- 
file derived from an NFW density profile to profiles mea- 
sured for isolated galaxies in the mock catalogs with the 
same redshift and magnitude range and search criteria as 
the data. The mock catalogs employed here have total 
volume that is comparable to the final DEEP2 sample. 
The solid line is the profile from the mock catalogs be- 
fore we include the effects of target selection. The best- 
fit NFW profile (dotted line) is derived as described in 
§ 13 The halo mass associated with this profile (M200) 
is 3.8io g7 X IO^^H-^Mq. We can derive the true halo 
mass for these host galaxies by using the simulations, 
since we know the number of dark matter particles in 
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Fig. 6. — The best fit NFW line-of-sight velocity dispersion pro- 
file (dotted line) to the derived profile for 480 satellites from mock 
catalogs (solid with error bars). The mass associated with this 
NFW profile is 3.8 X 10^^h~^ Mq, which agrees quite well with the 
true average halo mass of these galaxies, 4.2 X 10^^ Mq. The 
dashed profile with error bars was derived from the mock catalogs 
after including target selection effects. The errors increase due 
to a smaller sample size and contamination from galaxies falsely 
identified as isolated. The best-fit NFW curve for this sample cor- 
responds to a mass of 3.9 X 10^^ Mq, and again accurately re- 
covers the true halo mass. Note the different y-axis scales between 
this and Figs. 0and|5] 



the halo and hence can directly compute Af2oo- The true 
halo mass distribution for isolated galaxies in the mock 
catalogs conforms to a log-normal distribution with la 
range 3.0 x lO^^h'^MQ to 1.0 x IO^/i-^Mq with a peak 
of 4.2 X 10^^h~^MQ. The velocity profile fit therefore 
recovers the true average mass quite accurately, giving 
us confidence that the algorithm used to fit the profile 
robustly recovers the underlying halo mass. 

The dashed line in Fig. shows the velocity disper- 
sion profile derived from mock catalogs to which the 
DEEP2 target selection algorithm (which will obtain 
spectroscopy of only ~ 65% of eligible objects) and 70% 
redshift completeness have been applied. This profile is 
noisier both because of mistaken identification of galaxies 
as isolated (due to their neighbors not being included in 
the sample) and a 50% decrease in the number of satel- 
lites due to the combined effects of target selection and 
redshift incompleteness. 

The mass associated with the best-fit NFW profile is 
3.9li J X 10^2/i"^Mq. The true masses of the host ha- 
los conform to a log-normal distribution with la range 
1.7 X lOi2/i-iM0 to 5.7 X lO^^h-^MQ and peak at 
2.6 X 10^^h~^ Mq. When we compare the isolated hosts 
found in the pre-target selection mock catalogs to those 
found after target selection, we find that ^ 30% of hosts 
in the mock catalogs after target selection are not truly 
isolated, but only appear isolated because their compan- 
ions were not targeted for observation. Fortunately, this 
level of contamination does not seem to limit our ability 
to accurately recover the underlying halo mass, as to first 
order it mimics the effects of background interlopers and 
hence is accounted for by our interloper correction. 
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6. DISCUSSION 

We find that isolated galaxies at z ^ 1 have a simi- 
lar halo mass as isolated galaxies which are 1 magnitude 
fainter at z ~ 0. This implies that there has been little 
or no evolution in the halo mass of isolated galaxies with 
magnitudes in the range ~ M^-0.5 to ~ even 
though M ^ has evolved by ^ 1 magnitude over this red- 
shift range. Our results are thus consistent with isolated 
galaxies of a fixed luminosity relative to M* populat- 
ing their dark matter halos in a similar way from z ^ 1 
to z ^ 0, a result at tained with an independent method 
bv lYan. Madewick. fc White ( 2003). Phrased differently, 
we find that the dynamical _B-band mass-to-light ra- 
tio {M2oo/Lb) is increasing from M200/LB = 104 ± 43 
/iMq/Lq,s at z ~ 1 to M200/LB = 257±52 HMq/Lq^b 
at z ~ 0, a factor of 2.5. This increase is attributable 
solely to the 1 magnitude decrease in the typical satellite 
host galaxy luminosity from z ~ 1 to z ~ 0. Assuming 
that the isolated galaxies found in DEEP2 at z 1 pas- 
sively evolve to the SDSS isolated galaxies at z ~ 0, our 
results imply that the ratio of baryonic mass to dark halo 
mass in these galaxies has been constant for the last '--^8 
billion years. 

These results are insensitive to the various assump- 
tions used to calculate masses; when we instead assume 
an isothermal model for the dark matter density distri- 
bution, the values of the masses calculated for both the 
high and low redshift samples change, but the consis- 
tency between these two values remains. Specifically, for 
the isothermal model our isolated galaxy halo mass at 
z 1 within r2oo becomes 1.9 ± 0.95 x 10^^ H'^Mq 
while for isolated galaxies at z ~ th e mass within r^nn 
changes to 2.0±0.56xl0^^ H'^Mq fsee lBrvan fc NormanI 
[199 8. for the relevant relation between velocity disper- 
sion and mass within r2oo in an isothermal model). The 
same can be said for our other assumptions, including 
the anisotropy required in the Jeans equation and the 
concentration for the NFW profile. The value of the con- 
centration parameter c is actu ally expected to be lower 
at high redshift than z ~ l|Bullock et al.U200i]|) . but 
this, too, should not have a strong effect. For instance, 
if c = 5 at z '--^ 1 instead of 10 f following IBuUock et al.l 
120011 who find that c cx (1 -I- z)~^ for halos of the same 
mass), M200 increases by ~ 1 x 10^^ h~^MQ, well within 
the quoted errors. It is essential for this comparison that 
masses at low and high redshift be calculated by the same 
method and with the same assumptions, but the specific 
assumptions do not affect our overall conclusions. 

Host-satellite kinematics in the local universe have 
also been studi ed in the 2dF Galaxy Redshift Survey. 
iBrainerdI l)2005tl measure a velocity dispersion profile for 
2L* isolated galaxies and find that it falls from ^ 200 
km s-i at lOO/i-i kpc to ~ 160 km s'^ at AOOh'^ kpc, 
in good agreement with what we measure here for SDSS 
galaxies of similar luminosity. 

Importantly, the halo mass derived from the data 
agrees with the mass derived from the mock catalogs for 
the same host magnitude range, —22 < Mb — 5 log(/i) < 
—21. Requiring such agreement can be used to set con- 
straints on the way in which the number of galaxies in 
a dark matter halo and their luminosities depend on the 
underlying halo mass - i.e., the halo occupation distri- 
bution and conditional luminosity function, which were 



used to create these mock catalogs. 

Comparison to other methods for determining halo 
masses, such as galaxy-galaxy lensing, is complicated for 
a number of reasons. First, to date galaxy-galaxy lens- 
ing studies have focused on ^ L* galaxies, while here 
we study a population of galaxies with mean luminosity 
~ 2.5L*. Second, galaxy-galaxy lensing probes all of the 
mass along the line-of-sight, not just the mass around an 
isolated galaxy. Third, we require satellites around iso- 
lated galaxies in order to determine the halo mass while 
galaxy-galaxy lensing can probe galaxies without satel- 
lites. Halos with luminous satellites could potentially 
be systematically more massive than halos without lu- 
minous satellites. Yet even with these potentially im- 
portant differences, the evolutionary trend described in 
this paper is in good agreement with galaxy-galaxy lens- 
ing results which show that the halo mass of isolated 
~ L* galaxies is approximately constant from z ^ 0.8 to 
z ^ .15 (Wil son et al.ll200lt iHoelcstra. Yee. fc Gladdersl 
l200llKleinheinrich et alJl2005j) . 

As mentioned in § 0] our results at z ~ are in good 
agreement with P03 when we use their definitions of host 
galaxies and masses. This good agreement is encourag- 
ing and implies that our mass reconstruction method is 
robust, as we do not exactly follow the mass estimation 
method outlined in P03. Mass estimates at z ~ for 
isolated galaxies utilizing satellite kinematics seem to be 
converging. 

It is important to keep in mind that, due to the small 
number of satellites found, this analysis has been applied 
to a host sample consisting of both early and late type 
galaxies, and is therefore mixing two different popula- 
tions of galaxies. We thus stress that these are initial 
results requiring more data to untangle these sorts of 
complications. 

Upon completion of the DEEP2 survey we will have a 
sample 4x larger than what was used for the present 
analysis, which will decrease our uncertainties on velocity 
dispersions by a factor of 2. This will result in a similar 
decrease in errors on our halo mass estimate, allowing 
for much tighter constraints on both halo mass evolution 
and the halo model. With the completed survey, we will 
be able to separate our host galaxies by spectral type, 
color, or redshift, allowing for more precise comparisons 
to local samples. 
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TABLE 1 

SEARCH PARAMETERS AND DERIVED QUANTITIES FOR 7 SAMPLES. S REFERS TO 
SATELLITE SELECTION CRITERIA, WHILE A REFERS TO THE PARAMETERS FOR 
HOST GALAXY ISOLATION. 



Sample SMb Srp Sv AMb Arp Av N ^ V200 M " 



1 


1.5 


350 


500 


0.75 


700 


1000 


75 


41 n+^^ 


5.5 


2 


1.0 


500 


750 


0.75 


1000 


1000 


160 


4101^5 


5.4; 


3 


1.5 


500 


750 


1.0 


1000 


1000 


52 


5.3 


4 


1.5 


350 


700 


1.0 


700 


1000 


55 


410l- 


5.5 


5 


1.0 


500 


700 


1.0 


1000 


1000 


11 


400l™ 


5.0: 


6 


1.5 


500 


500 


1.0 


1000 


1000 


51 


4251^7^ 


6.1; 


7 


1.5 


350 


500 


1.0 


500 


1000 


82 




6.1; 



+2.5 



1+3.1 



''Total number of satellites with host magnitude —22 < Mb — 5log{h) < —21. 

''V200 measured in km 

'^Maoo/lOi^ft-iMQ 
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